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GSFC & JPL partnership 

detect low frequency GW 
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LIGO: Hanford, WA, and 
Livingston, LA; L = 4 km 

VIRGO: PISA, L = 3 km 
GEO600: Hannover L = 600 m 
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Complementary observations, different frequencies & sources.. 


► compact systems 
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Simulation of th 
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Ground-based detectors 

BH binary coalescences 
#■/ 



direct 


LISA can confront GR with (nasa/cxc/mpe/. 

observations in the dynamical, strong field regime . . . 

.... if we know the merger waveforms 






GWs from final merger of Mack hole Unary... 


lute 


to calculate dynamics & waveform! 
spins ^ apply to ground-based & I 


(graphic courtesy of Kip Thorne) 
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a major theory 
challenge 
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Nearly as difficult as building these (gravitational wave j 
observatories, however, is the task of computing the 
gravitational waveforms that are expected when two black holes 
merge. This is a major challenge in computational general 
relativity and one that will stretch computational hardware and 
software to the limits. However, a bonus is that the waveforms 
will be quite unique to general relativity, and if they are 
reproduced observational ly, scientists will have 
performed a highly sensitive test of gravity in 
the strong-field regime.” 




Constraint equations 
Evolution equations 


lapse function shift vecto 
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A Brief History of BBH simulations. 




1970s: Smarr an J 171 


“wormholes” 
ppley: head-on collision of 2 BHs, extract GWs 


Pioneering efforts on supercomputers at Livermore Natl Lab 


I M 


Work on 2-D head-on collisions at NCSA 

NSF Grand Challenge: multi-institution, multi-year effort in 3-D 


^ This is reallv 


'fficult! Instabilities 9 issues in formalisms, etc. 


Diaspora: multiple efforts (AEI, UT-Austin, PSU, Cornell...) 
Difficulties proliferate, instabilities arise, codes crash.... 


xerical relativity is 9 
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Os: LIGO/GEO/VIRGO and LISA spur more development 

M. JL 

New groups arise: Caltech, UT -Brownsville, LSU, NASA/GSFC. 


Breakthroughs & rapid progress throughout community 
Orbits, waveforms, and astrophysical applications.... 



Issues and ingredients for success... 

Formulations of the Einstein equations 

- fully 2 nd order, fully 1 st order, mixed 1 st and 2 nd order PDEs 

- which variables to use? 

- incorporate constraints into evolsi eqns? solve constraints? 

Coordinate conditions: 

- lapse function __ - “singularity avoiding” time slicing 

- shift vector - keep coordinates from falling into black holes., 

Constrained initial data to approximate astrophysical binary 

- start on approx quasi-circular orbits 

- inward radial velocity... 


w to handle the black holes: 

excision? punctures? 

comoving coordinates? move the black holes? 


ananiA crrirf rpcsi 


sn ci • 


~ (10 - 100)M 

c = G = 1^1M~5x 10" 6 (M/M Sun ) sec - 1.5 (M/M Sun ) km 
finite differences w/ mesh refinement; spectral methods 



Apparent Hbrizon 1 and Lapse at Time T~ 1 14M 


orbital plane 


^ fix liie BH iimicliires in the grid 


Cl i 


> Tich y, & Jansen, PR] 
2004), gr-qc/G3l21 12 

nonspinning BHs 
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The 1 st orbit, merger, & ringdown 


- metric g V} is basic variable 

- 2 nd order in space & time 

Incised BHs move through gj 
AMR: high resolution arounc 


edge of g 

o O 




M - E, l. v -40M! 


new idea: moving puncture BHs 


Campanelli, et al., PRL, 96, 
111101 (2006), gr-qc/0511048 

Baker, et al., PRL, 96, 111102 
(2006), gr-qc/0511103 


Regularize near puncture 
New conditions for & ' 
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A powerful new idea 


Conformal formalism, BHs represented as punctures 


Van Meter, et al., “How to move a puncture black hole without 
excision...,” PRD 73 (2006) 124011 2006), gr-qc/0605030 


0.004 


Last Orbit 

I SCO (transited b# 95M) 


0.002 


PSU, AEI/LSU, FAU/Jena... 

At April 2006 APS meeting, a 
full session devoted to BBH 
mergers w/ moving punctures! 

Summer 2006: method adopted 
by most of community 


0.002 


0.004 


Campanelli, etal., PRD, 73, 061501 
(2006), gr-qc/06010901 
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Campanelli, et al., gr-qc/0612076 




<er, ct ah, ApJL, 653, L93 (2006) astro-j3.il/0603204: 


q = m 1 /m 2 = 0.67 

widest separation run completes 
- 2.5 orbits before merger 

agrees w/ PN over most of 
1 st orbit to better than 1% 

Overall, report kick values in the 
range v kick = (86 - 97)km/s 

mzales, ct al, gr-qc/06 1 0 1 54 
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Ran series of runs, w/ mass ratios 
in the range 0.253 < q < 1 
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Recoiling from mergers of spinning BHs 


Herrmann, et aL, gr-qc/0701 143 


Koppitz, et aL, gr-qc/0701 163 


f « W 


Campanelli, et al., gr-qc/0701 164 


Gonzalez, et al., gr-qc/0702052 


Baker, et al., astro-ph/0702390 


interesting parameter space... more studies to come 



Current status of BBH merger simulations... 

Impressive recent progress on a broad front: many research groups, 
different codes, methods. . . 

Equal mass, nonspinning BBHs: several groups are now evolving 
for several orbits, followed by the plunge, merger, and ringdown 

There is general agreement on the simple waveform shape and that 

- total GW energy emitted in last few cycles 

_E ~ (0.035 - 0.04)M (depends on the number of orbits) 

~ final BH has spin a ~ 0.7M 

Long runs now possible. . 7 orbits before merger 

Applications to G W data analysis are beginning 

Explosion of work on nonequal mass and spinning BH mergers 
and the resulting kicks 

- Interesting parameter space 

- Important astrophysical applications. . . 
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simulations 



